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ABSTRACT 


A  simple  method  for  the  rapid  determination  of  the 
linearized  surface  pressure  distribution  on  streamlined  ducted 
bodies  in  axisymmetric  flow  is  developed  using  conventional 
singularity  distributions.  The  final  form  of  the  pressure 
coefficient  is  expressed  directly  as  a  single  matrix  operation 
on  the  duct  sectional  properties  and  tables  of  the  necessary 
matrix  coefficients  are  given.  Several  numerical  examples  are 
worked  out  and  compared  with  the  results  of  previous  theoretical 
and  experimental  investigations. 
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THE  AERODYNAMIC  LOADING 
ON 

STREAMLINED  DUCTED  BODIES 


INTRODUCTION 


During  the  past  three  years  we  have  formulated1  and 
studied  ' J  a  three-dimensional  theory  for  the  f inite-bladed 
shrouded  propeller  in  the  forward  flight  regime.  The  mathematical 
model  chosen  to  represent  the  system  was  based  on  the  classical 
concepts  of  vortex  propeller  theory  and  of  lifting-surface  theory 
for  ring  wings  of  zero  thickness  in  inviscid,  incompressible  flow. 
The  purpose  of  this  report  is  to  remove  the  limitation  of  zero 
thickness  as  all  shrouds  have  a  finite  thickness  which  may  con¬ 
tribute  substantially  to  the  determination  of  such  quantities  as 
the  shroud  loading  and  propeller  inflow,  since  the  effect  of  the 
propeller  may  be  simply  superimposed  within  the  framework  of 
linearized  theory,  it  is  sufficient  for  our  purposes  to  limit  our 
attention  directly  to  the  problem  of  the  annular  airfoil. 

The  problem  of  finding  the  detailed  velocity  field  associated 
with  and  the  forces  acting  on  annular  airfoils  in  incompressible 
flow  has  been  the  subject  of  extensive  investigation  over  the 
past  twenty-five  years,  possibly  the  earliest  treatment  using 
singularity  distributions  and  classical  thin-airfoil  theory  was 

given  by  H.  E.  Dickmann^  in  1940.  Since  that  time,  many  contri- 

R-12  1^-1Q 

butions^  have  appeared,  culminating  in  the  studies 
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of  J.  Weissinger,  Who  was  the  first  to  give  a  comprehensive 
treatment  of  the  thickness  problem.  Nevertheless,  in  spite  of 
all  the  work  which  has  been  undertaken,  a  simple,  accurate  method 
had  not  yet  been  developed  for  the  rapid  prediction  of  duct 
surface  pressure  distributions  for  preliminary  design  studies. 

In  this  report,  such  a  method  for  the  linearized  pressure 
distribution  in  the  case  of  axisymmetric  flow  is  developed.  The 
results  obtained  are  presented  in  a  form  such  that  they  may  be 
incorporated  directly  into  our  previous  studies  of  the  ducted 
propeller  or  used  "as  is"  in  the  treatment  of  ducted  bodies. 

The  conventional  mathematical  model  of  ring  vortices  and 
sources  distributed  on  a  reference  cylinder  is  chosen  to  repre¬ 
sent  the  duct  and  the  approximations  of  thin-airfoil  theory  are 
employed  throughout.  The  source  strength  depends  only  upon  the 
rate  of  change  of  the  airfoil  profile  thickness.  In  turn,  this 
distribution  induces  a  radial  wash  which,  together  with  that 
induced  by  the  vortex  distribution,  must  be  equal  to  the  local 
camber  everywhere  on  the  duct  surface,  with  the  thickness  dis¬ 
tribution  given  in  a  generalized  form  of  the  NACA  four -and  five¬ 
digit  airfoil  series,  the  thickness -induced  radial  wash  is 
expressed  in  appropriate  matrix  form.  The  solution  for  the 

unknown  vortex  strength  follows  immediately  using  the  iteration 

o 

technique  of  D.  E.  Ordway  and  M.  D.  Greenberg  .  Then  the 
linearized  surface  pressure  distribution  is  reduced  to  a  single 
matrix  operation  on  the  duct  sectional  properties.  The  effect 
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of  the  duct  chord-to-diameter  ratio  is  explicitly  obtained  and 

tables  of  all  necessary  matrix  coefficients  are  given  for  four 

typical  values  of  this  parameter.  Incidentally,  several  defi- 

20  21 

nite  integrals  not  appearing  in  common  integral  tables  ' 
were  evaluated  during  the  course  of  this  work  and  have  been 
collected  in  an  appendix. 

Briefly,  the  report  is  presented  as  follows.  The  first 
Chapter  gives  the  mathematical  formulation  of  the  problem,  the 
derivation  of  the  induced  velocity  field,  and  the  procedure  for 
finding  the  strengths  of  the  source  and  vortex  distributions . 
Chapter  Two  details  the  reduction  of  the  thickness -induced 
radial  wash  to  matrix  form  in  terms  of  the  known  thickness  dis¬ 
tribution  coefficients.  In  Chapter  Three,  the  duct  surface 
coefficient  is  found  in  matrix  form  and  simple  expressions  for 
the  sectional  radial  force  and  pitching  moment  coefficients  are 
furnished.  Tabulated  values  of  all  the  required  matrix  coeffi¬ 
cients  for  four  body  chord-to-diameter  ratios  are  presented. 
Several  numerical  examples  are  worked  out  to  illustrate  the 
computational  technique  and  to  compare  with  previous  theoret¬ 
ical  and  experimental  results  in  Chapter  Four. 


CHAPTER  ONE 


BASIC  FORMULATION 


1.1  Mathematical  Model 

Consider  an  annular  airfoil,  or  duct,  at  zero  incidence  in 
a  uniform,  inviscid,  incompressible  stream  u  ,  see  Fig.  1.1. 

A  cylindrical  coordinate  system  (x,r ,0)  is  chosen  with  the 
origin  fixed  at  the  center  of  the  duct,  x  =  -c/2  corresponding 
to  the  duct  leading  edge  and  x  =  +c/2  to  the  trailing  edge. 

The  duct  section  is  assumed  known  and  is  the  same  at  every 
azimuthal  station.  The  local  inclination  of  the  camber  line  to 
the  x-axis  is  denoted  by  e(x)  and  the  thickness  distribution 
by  t(x)  .  Both  e  and  t/c  are  assumed  small. 

Making  the  conventional  approximations  of  linearized  per¬ 
turbations,  the  annular  airfoil  is  replaced  by  a  distribution  of 
ring  sources  and  vortices  on  an  appropriate  reference  cylinder 
of  radius  R  .  Due  to  the  axisymmetry  of  the  flow  field,  the 
source  strength  per  unit  length  f  does  not  vary  with  6  . 

The  same  is  true  of  the  vortex  strength  per  unit  length  y  . 
Since  7  =  >(x)  ,  there  are  no  trailing  vortices  from  the  duct. 

The  strength  of  these  distributions  will  be  determined 
from  the  physical  boundary  condition  that  the  flow  be  everywhere 
tangent  to  the  duct  surface  together  with  the  Kutta-Joukowski 
condition  at  the  trailing  edge,  consistent  with  our  assumptions, 


4 


6 


this  reduces  to  the  familiar  thin-airfoil  approximation 


t  = 


vf  +  vy 

u 


(1.1) 


where  vf  and  v^  are  the  radial  velocities  induced  on  the 
reference  cylinder  by  the  source  and  vortex  distributions 
respectively  and  7 (c/2)  =  0  .  Throughout  the  report  we  adopt 
the  ccnvention  that  the  prime  (')  denotes  total  differentiation 
of  a  function  with  respect  to  its  indicated  argument  and,  in 
the  event  a  double  sign  precedes  a  quantity,  the  top  sign  is  to 
be  used  to  give  the  value  on  the  outer  duct  surface  and  the 
bottom  sign,  the  value  on  the  inner  surface. 


1.2  induced  velocity  Field 

To  solve  Eq.  (l.l)  the  radial  component  of  the  velocities 
induced  by  f  and  7  are  required.  Subsequently,  the  axial 
component  will  be  needed  for  the  pressure  coefficient.  There¬ 
fore,  we  will  derive  general  expressions  for  both  before  we 
proceed . 

The  velocity  potential  <t>  for  an  isolated  source  ring 
of  unit  strength  is  generally  expressed  in  terms  of  the  complete 

elliptic  integral  of  the  first  kind.  Alternatively,  we  may  use 

22 

the  relation  between  this  function  and  the  Legendre  function 
of  second  kind  and  minus  one-half  order  to  obtain, 
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«X>  -  -  -i—  q  ,  (m) 

o>  s  i  +  ^±i£=ili  (1.2) 

2r 

where  x/R  =  x  ,  r/R  =  r  ,  ...  and  Ax  li  the  axial  separation 
of  the  field  point  and  the  singularity  at  i  ,  Ax  =  x  -  £  • 

The  corresponding  axial  and  radial  velocities  are  gotten  from 
d<$/dx  and  b<t>/bz  respectively.  Upon  reduction  of  the  latter 

pp 

velocity  by  recursion  relations  and  chordwise  integration  of 
both  over  f  ,  we  find  for  the  resultant  contribution  uf  to  the 
axial  perturbation  velocity, 

X 

uf(x,r)  -  — J  f (C)  Ax  Q'^(oo)  d£  (1.3) 

27rr  -X 

and  v£  to  the  radial  perturbation  velocity, 

X 

vf(x,r)  »  — f  fQ^(®)  "  d$  U-2*) 

“  X 

where  X  =  c/2R  and  Q^(cb)  is  the  Legendre  function  of  second 
kind  and  plus  one-half  order. 

On  the  other  hand,  the  velocity  components  for  an  isolated 
vortex  ring  of  unit  strength  may  be  found  most  easily  from  the 
stream  function  Y  .  in  terms  of  Q^(a>)  ,  we  have^  t 

*  -  ^  ^  V®*  (l,5) 
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X 

v7(x,l)  m  f  7(4)  Ax  Q^(Z)  d£ 

-X 

where 


(1.11) 


(1.12) 


1.3  Determination  of  Source  and  vortex  Strengths 

Formally,  f  and  y  are  coupled  by  Eq.  (1.1).  However, 
from  Eqs.  (1.9)  and  (l.ll),  only  vf  and  not  is  discontin¬ 

uous  across  the  reference  cylinder.  Thus  we  must  have 


f(x)  =  Ut'(x) 


(1.13) 


and  the  required  source  strength  depends  only  upon  the  rate  of 
change  of  the  duct  thickness  distribution  as  for  the  problem  of 
the  two-dimensional  airfoil10. 

Consequently,  with  f  determined  we  may  regard  the  contin¬ 
uous  part  of  v^  as  a  thickness -induced  camber  defined  by 

X 

«t  s£rf  (Q^(S)  -  Q;ij(S))  4  (1.14) 

-  X 

and  combined  with  €  to  form  an  effective  camber  €e  , 


(1.15) 


10 


Then  Eq.  (1.1 )  can  be  rewritten  as 

X 

€e  "  "  f  A*  o£(®)  (1.16) 

-X 

where 

7  H  UC(x)  (1.17) 

Eq.  (1.16)  constitutes  an  integral  equation  for  the  unknown 
vortex  distribution  C  .  To  solve,  we  shall  use  the  iterative 
procedure  of  Ordway  and  Greenberg.  A  new  axial  coordinate  0 
is  introduced  such  that 

x  =  -  X  cos  0  (1.18) 

where  0  <  0  <_ n  and  c  is  expressed  in  a  Glauert  series, 

00 

C  =  cQ  cot  §■  +  X  cv  sin  v<t>  (1.19) 

V  =1 

Such  a  series  satisfies  the  Kutta-Joukowski  condition  at  the 
trailing  edge  and  possesses  the  proper  square-root  singularity 
at  the  leading  edge  in  the  cotangent  term.  The  corresponding 

p 

Glauert  coefficients  are  then  given  in  matrix  form  as 

(c)  *  (tlJ  +  [p]  +  [p]2  +  [p]3  +  ...  )(c2_D)  (1.20) 


where 
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(1.21) 


and  (c2-D)  is  a  column  matrix  of  the  Glauert  coefficients  of 
a  two-dimensional  airfoil  having  a  camber  equal  to  ,  these 
coefficients  being  given  by. 


c2-D.2  r  c  m 

0  it  J  e  r 


7T 

2—d  4  r 

Cy  ■  -  1  J  e  e  cos  V0  d(p  v  >  1 

0 


(1.22) 


['I-J  is  the  unit  matrix  and  [p]  is  a  square  matrix  whose 
coefficients  are  functions  only  of  X  .  Each  term  of  the  operator 
(flO  +  [p]  +  [p]  +  ...  )  corresponds  to  a  step  in  the  iteration 

process  which  in  the  limit  converges  to  the  exact  solution,  cf. 
Ref.  13.  The  leading  forty-nine  elements  of  [p]  are  given  in 
Ref.  2  for  X  =  1/4,  1/2,  3 /4  and  1  which  covers  the  range  of 
practical  interest.  For  convenience,  we  have  computed  the 

p 

elements  of  (['!-]  +  [p]  +  [p]  +  ...  )  for  these  X  to  an 
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_2i 

accuracy  of  +  1  x  10  .  These  are  presented  in  Tables  1.1- 

1.4.  Since  the  elements  diminish  rapidly  off  the  diagonal  and 
approach  unity  quickly  on  the  diagonal,  engineering  accuracy 
can  be  achieved  in  a  few  steps,  of  course  as  the  chord-to- 
diameter  ratio  increases,  the  three-dimensional  effect  increases 
and  the  speed  of  convergence  decreases. 

In  summary,  the  source  strength  is  found  immediately  from 
the  known  thickness  distribution,  see  Eq.  (1.13),  and  so  . 
From  €  and  e  ,  the  Glauert  coefficients  of  the  vortex  dis¬ 
tribution  are  then  found  from  a  single  matrix  operation  on  the 
corresponding  coefficients  of  Eqs.  (1.22)  for  a  two-dimensional 
airfoil  with  camber  € 

e 
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TABLE  1.4 

ELEMENTS  OP  (flj  +  [pj  +  [p]2  +  ...  )  FOR 


CHAPTER  TWO 


REDUCTION  OF  EFFECTIVE  CAMBER 


2.1  General  Considerations 

From  Eqs.  (1.22)  we  see  that  the  determination  of  the 
two-dimensional  Glauert  coefficients  is  equivalent  to  decom¬ 
posing  ee  into  the  Fourier  cosine  series 


e 

e 


cos  v$ 


(2.1) 


i  •  6  •  i 


(c2_D)  -  {c#]  5  C € )  -  (€fc)  (2.2) 

We  will  assume  that  (e)  or  the  contribution  of  the  geo¬ 
metrical  camber  to  (ee)  is  known  by  evaluation  of  Eqs.  (1.22) 
with  replaced  by  e  ,  cf.  Ref.  24,  or  by  direct  expression 

of  €  in  the  form  of  Eq.  (2.1),  cf.  Ref.  25 •  The  contribution 
of  the  thickness-induced  camber  will  now  be  determined. 

To  carry  this  out,  t'  and  the  second  factor  in  the  integrand 
of  Eq.  (1.14),  say  l^  ,  are  expressed  in  forms  suitable  for 
integration.  We  find  each  of  these  is  composed  of  a  singular 
as  well  as  a  regular  part  and  the  corresponding  products  are 
integrated  in  turn. 
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2.2  Form  of  Thickness  Distribution  and  t' 

We  will  assume  that  t  is  expanded  about  the  leading 
edge  in  a  generalized  form  of  the  NACA  four-  and  five-digit 
airfoil  series2^ 


00 


where  Aq  ,  A^  ,  Ag  . . .  are  known  constants  and  finite  in 
number  for  the  NACA  four-  and  five-digit  sections  and  may  be 
found  by  appropriate  curve  fitting  for  others .  The  leading 
coefficient  AQ  is  related  to  the  airfoil  leading -edge  radius 
1  by 


This  form  of  the  thickness  distribution  has  a  square-root 
behavior  at  the  leading  edge  only,  cf.  Refs.  11  and  15  in  which 
the  assumed  distribution  possesses  inherently  not  only  a  square- 
root  behavior  in  t  at  the  leading  edge  but  also  at  the  trail¬ 
ing  edge  as  well. 

Differentiating  Eq.  (2.3)  and  rewriting  in  a  power  series 
in  x  for  convenience,  we  find 
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Where  the  respective  coefficients  are  related  by 

a° -/K 

al=Al  +  A2+£A3+l‘A4  +  &A5+&A6  +  ^rA7+--- 

a2  “  r<A2  +  2A3  +  2A4  +  ?'A5+li‘A6+fi‘AT  +  *  *  *  ^ 

a3  “  a3  +  i  A4  +  F  a5  +  ¥*  a6  +  T?  *7  +  ‘  *  *  } 

a4  "  a4  +  I  a5  +  ¥  a6  +  l£  *7  +  *  *  •  ) 

a5  "  a5  +  ll  a6  +  ^  *7  *  *  *  } 

a6  *  ^5^a6  +  §5“  *7  +  ) 

®7  “  *7  +  *“  ) 

(2.6) 

For  NACA  four-  and  five-digit  airfoils,  these  coefficients  are 
tabulated  in  Table  2.1.  From  Eq.  (2.5)  then,  we  see  that  t' 
is  composed  of  a  square-root  singularity  at  the  leading  edge  and  a 
regular  part. 


0.20944  0.29690  0.36363  0.41988 
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2.3  Form  of  1^ 


Similarly,  the  second  factor  in  ,  or 

zi  s  b  {q4(£)  -  °V£))  (2*7) 


may  be  separated  into  two  parts,  a  logarithmic  singularity  and 

PP 

a  regular  term  (R^  .  From  the  known  behaviour  of  and 

Q  ^  near  unit  argument,  we  find 


-2 


In  +  Ol1(Ax) 


(2.8) 


and  both  the  singular  and  regular  parts  are  symmetric  about  Ax  *  0  . 
The  regular  part  has  been  evaluated  and  is  shown  in  Fig.  2.1. 

Since  (R^Ax)  =  (R^X  cos  cp  -  X  cos  <p)  is  bounded,  it  is 
convenient  to  expand  it  in  a  double  Fourier  cosine  series  in  the 
region  0  <_  0  <  7r  and  0  <  cp  <  tt  (C  =  -X  cos  cp)  as 


00  00 

&1  =  I  I  hp  q(X)  cos  cos  qcp 
p=0  q=0 

The  coefficients  h  (X )  are,  from  orthogonality, 

P  t  SL 


(2.9) 


n  tt 


*P,q  =  ^ 


SI*- 1 


cos  p0  cos  q <p  d<p  d</> 


(2.10) 


0  0 


for  p,q  /  0  ;  hp  Q  and  h0  q  are  one-half  and  hQ  Q  is 
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one-quarter  of  these  values.  Since  is  an  even  function 

in  Ax  ,  if  we  interchange  <f>  and  q>  in  Eq.  (2.10)  we  can 
show 


p.q 


=  h 


q,p 


(2.11) 


In  addition,  replacing  <t>  by  (ir-$)  and  «p  by  (w-<p)  in 
Eq.  (2.10),  we  find 


h  =  0  if  (p+q)  =  odd  integer  (2.12) 

P  t  vl 

For  X  =  1/4  ,  1/2  ,  3/4  and  1  we  have  computed  hp  q  for 
p,q  =0,  1,  2  ...7  from  Eq.  (2.10)  using  a  relaxation  proce¬ 
dure  developed  by  G.  E.  Bartholomew*^.  This  was  carried  out 
until  the  difference  between  successive  computations  was  less 
than  +  1  x  lO-'’  and  the  results  are  presented  in  Tables  2.2- 
2.5. 

A  comparison  of  the  exact  value  of  (R^  as  obtained  from 

28 

Eqs .  (2.7),  (2.8)  and  tables  of  the  Legendre  functions  with 

the  series  value  as  given  by  Eq.  (2.9)  using  the  computed  ^p^'8 

-4 

is  given  in  Table  2.6.  We  find  agreement  to  within  +  1  x  10 
in  all  cases,  the  results  becoming  significantly  better  as  X 
decreases .  This  is  due  to  the  rapid  decay  of  the  h^  ^  with 
increasing  p  and  q  at  low  values  of  X  .  As  expected,  the 
use  of  more  terms  in  the  double  Fourier  series  gives  a  higher 


COEFFICIENTS  h  „  FOR  X  =  3/4 


.04920  0.04921  0.04920 
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degree  of  accuracy. 

An  independent  check  on  the  accuracy  of  the  relaxation 
technique  was  desired.  Accordingly,  by  the  change  of  vari¬ 
ables  y  =  x  +  f,T}=x-£,we  found  it  possible  to  perform 
analytically  the  integration  on  y  and  so  reduce  h  to 
the  evaluation  of  a  single  integral.  Unfortunately,  except 
for  hQ  q  ,  the  complexity  of  the  integrand  precludes  the 
possibility  of  further  numerical  computation;  hQ  q  becomes 

o  P 

h0, 0  -  ^  J  -W=T  Q-h  {1  +  ^(2X-t,)-2)  dT,  (2.13) 

This  was  evaluated  for  X  =  1/2  and  1  and  the  results  agree 
to  within  one  unit  in  the  fifth  decimal  place  with  those 
determined  by  the  relaxation  method. 

With  the  expansion  of  both  t'  and  into  a  singular 

and  a  regular  part,  we  may  proceed  to  the  determination  of 
(et)  .  The  integral  form  of  ,  Eq.  (1.14),  can  be  sepa¬ 
rated  into  four  components  corresponding  to  the  multiplication 
of  the  singular -singular ,  singular -regular,  regular-singular, 
and  regular-regular  parts  of  t'  and  1^  respectively.  We 
will  now  consider  the  evaluation  of  each  of  these  in  turn  and 
the  subsequent  decomposition  of  (et)  into  the  matrix  form 
[q]  (a)  . 
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2.4  Singular-Singular  Contribution 

The  first  contribution  to  arises  from  the  product  of 

the  two  singular  terms,  upon  integration  by  parts  and  trans¬ 
formation  to  the  Glauert  angular  variable  0  ,  we  find 


2  +  ~  fn( l+cos0 )  -  sin*s0  aQ  (2.14) 


The  quantity  inside  the  brackets  is  finite  for  0  <_  0  <  tt  as 
cam  be  verified  by  passing  to  the  appropriate  limits. 

To  obtain  the  contribution  to  (et)  ,  Eq.  (2.14)  is  sub¬ 
stituted  into  Eqs.  (1.22)  and  the  necessary  integrations 
carried  out  using  elementary  integrals  and 


in(l+cos0)  d 0  =  -tt  in  2 


cos  n0  in(l+cos0)  d0  = 


n 


sin  0  in 


l+sin0 

l-sin0 


d0  =  27r 


n  >  1 


(2.15) 


(2.16) 


(2.17) 


and 
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cos 2n<fc  sin  <(>  d0  = 


2n  1  •  3  •  5  »  •  •  (2n-l) 

2n+l  2  •  4  •  6  •  •  •  (2n) 


n  >  1 


(2.18) 


Eq.  (2.15)  is  taken  from  Ref.  20.  Derivations  of  Eqs.  (2.16), 
(2.17)  and  (2.18)  may  be  found  in  the  Appendix. 

With  these  results,  the  contribution  to  t3  may  be 
written  in  the  matrix  form  [Qgs](a)  where  [Qsg]  is  given 
in  Table  2 .7 • 


2.5  Singular -Regular  Contribution 

in  view  of  the  form  of  the  expansion  of  in  Eq.  (2.9), 

the  contribution  to  of  the  product  of  the  singular  part  of 

t'  and  the  regular  part  of  1^  reduces  immediately  to  a  Fourier 
cosine  series,  or 

X 

/ 

-X 


00  00 

fl  i 


{*+(  p= 


h  „  cos  p 0  cos  q®  = 

Pig 


p=C  q=0 


00 

I 

p=0 


cos  p0 


2  Vix 


00 

I 

q=0 


p,q 


(-D 


q+i 


/ 


(4q2-l) 


(2.19) 


To  carry  this  out,  we  have  changed  to  the  Glauert  variable  and 
used  the  relationship  sin  5j<p  =  •>/ %{ 1-cos  q>)  .  Thus  the  contri¬ 
bution  to  ( € ^}  may  be  written  simply  by  comparison  with  Eq.  (2.1). 


o  o  o 


o  o  o 


o  o  o 


o  o  o 


o  o  o 


o  o  o 


h|oo 


Hl=t 
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Similarly  to  the  preceding  contribution,  this  may  be  ex¬ 
pressed  as  [ Qsr3  C a)  with  [QgR]  given  in  Table  2.8. 


2.6  Regular -Singular  Contribution 


The  product  of  the  regular  part  of  t'  and  the  singular 
part  of  l1  combine  to  give 


00 

I 

n=l 


1 

87rn 


X  _n 


-2  f  df  +  Xn  In . -  (-X)n  In— - 

J  Ax  (X-x)2  (X+x)2 


-X 


an  (2.20) 


The  integrals  on  the  right-hand  side  have  been  evaluated 
analytically2^,  e.  g., 


=  x  in  -  2X 

X-x 


=  £  ln  _  2Xx 
X-x 


(2.21) 

The  reduction  of  Eq.  (2.20)  to  the  matrix  form  [QRg](a)  is 
quite  lengthy.  To  illustrate  we  will  consider  only  a  typical  case. 
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say  n  =  1  .  From  the  first  of  Eqs.  (2.21)  and  Eq.  (2.20), 
this  term  becomes,  after  changing  to  the  Glauert  variable  and 
rearranging, 


35 


X 

55F 


o 

2  +  in-w-  -  (1-COS0)  Jn(l-cos0)  -  (1+cos 0)  in(l+cos0)  a. 
\  1 


2 

If  we  replace  in(l-cos0)  with  in  2sin  J50  and  in(l+cos0) 
2 

with  in  2cos  *s0  and  expand  the  resulting  forms  by  means  of 
formulas  603.2  and  603«^  of  Ref.  20  we  get 


_X_ 

2tt 


^  cos  20 


-  cos  40  - 


I5!> 


cos  60 


al 


Again,  the  contribution  of  this  term  to  (€fc)  is  found  by  com¬ 
parison  with  Eq.  (2.1).  In  a  like  manner,  the  contributions  of 
the  other  terms  are  obtained  and  so  [qrs]  *  see  Table  2.9- 


2.7  Regular-Regular  Contribution 

The  final  contribution  to  {ct}  from  the  regular  parts  of 

t'  and  I1  is  determined  in  the  same  fashion  as  the  singular- 

regular  term.  Here,  this  requires  the  evaluation  of  integrals 

_n-l 

whose  integrand  is  of  the  form  £  cos  q<p  .  This  is  carried 
out  by  changing  to  the  Glauert  variable  and  replacing  cos11  -1* 
by  its  alternate  cosine  series,  cf.  formulas  404,  Ref.  20.  The 
resulting  integrals  follow  in  turn  from 
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Ut)  =  [Q]  Ca)  (2.24) 

The  leading  forty-nine  elements  of  [q]  have  been  computed  to 
an  accuracy  of  +  1  x  10-^  for  \  =  1/4  ,  1/2  ,  3/4  and  1 
and  these  are  presented  in  Tables  2.11-2.14.  Tables  2. 2-2. 5 
were  used  to  compute  the  elements  of  [QgR]  and  [Qrr]  ' 

In  summary  then,  we  have  from  Eq.  (2.2) 

{C2-D}  =  UJ  -  [Q]  (a)  (2.25) 

where  the  elements  of  (a)  are  related  to  the  coefficients  of 
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LO|t-'  U)|l-'  LO|f-  00|l— 


X  13  XI  \  2 

ho,o  +  15  ho,2  "  1C5  ho,4  "  315"  ho, 6  +  x 


1  13  XI  2 

h2,0  "  l3h2,2  +  l05  h2,4  +  315  h2,6  +  •••)  x 


*4,0  -  A-h4,2  +  Sh4,4  +  Iifeh4(6  +  •••)  xC 


(Fhl,l  -  ihl,3  _  ^hl,5  ‘  ^hl,7  +  •••)  X' 


(^h3,l  -  ih3,3  _  ^h3,5  ~  ^h3,7  +  •••)  X' 


(fh5.1  '  35  h5,3  "  ^h5,5  '  ^  h5,7  +  x' 


h6.0  *  T5h6,2  +T&h6A  +mh6,6  +  X‘ 


TABLE  2.10 
MATRIX  [Q  ] 
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the  thickness  distribution  by  Eqs.  (2.6).  consequently,  the 
contribution  of  shroud  thickness  has  been  simplified  to  a 
single  matrix  operation  on  these  related  coefficients. 


CHAPTER  THREE 


DUCT  PRESSURE  DISTRIBUTION, 
SECTIONAL  RADIAL  FORCE  AND  PITCHING  MOMENT 


Where  p  is  the  local  static  pressure,  p^  the  free  stream 
static  pressure  and  p  the  fluid  density.  Consistent  with  our 
assumption  of  small  perturbations,  we  use  the  linearized  form  of 
Bernoulli's  equation,  or 


P  -  P®  “  -  PU  (uf(x,l)  +  yx.l)}  (3.2) 

With  this  result  and  Eqs.  (1.8),  (1.10),  (1.13)  and  (1.17)* 

Eq.  (3*1)  reduces  to 

X 

cp  -  +  C(x)  +  J  t'(i)  Ax  0^(5)  d£ 

-X 

X 

-  k  f  c(«>  (Q^®)  -  QlisC®)}  *1  (3.3) 

-x 

and  the  net  duct  loading,  or  the  inner  surface  pressure  coefficient 


44 


45 


minus  the  outer  surface  pressure  coefficient,  is  simply 


cp3  -  -  2  C(x) 


(3.4) 


We  see  that  c  is  composed  of  a  discontinuous  part  pro- 

Jr 

portional  to  the  local  vortex  strength  and  a  continuous  part 
composed  of  a  term  due  directly  to  the  duct  thickness  and  a  term 
due  to  the  vortex  distribution,  or  equivalently,  the  effective 
camber.  Analogously  to  the  thickness-induced  camber,  the  con¬ 
tributions  to  the  continuous  part  may  be  reduced  to  a  matrix 
operation  on  the  duct  sectional  properties .  This  is  readily  seen 
for  the  second  integral  if  we  compare  the  integrands  of  v^(x,l) 
and  Uy(x,l)  ,  Eqs.  (1.9)  and  (1.10),  and  note  the  form  of  c(x)  , 
Eq.  (1.19).  For  the  first  integral,  however,  it  is  necessary  to 
put  the  second  factor  in  the  integrand,  say  l2  ,  in  a  suitable 
form. 


3.2  Form  of  I, 


We  now  consider  ig  ,  or 


z2  S  S'  QVS) 


(3.5) 


Physically,  Ic,  is  equal  to  the  axial  velocity  induced  on  the 
reference  cylinder  by  a  ring  source  of  unit  strength  having  the 
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same  radius.  From  the  known  logarithmic  behaviour  of  Q_jj 
near  unit  argument,  we  find  that  Ip  may  be  separated  into  two 
parts;  a  Cauchy  singular  term  and  a  regular  term  Gig  ,  or 

I  (dx)  -  -  -L_+  ft2(Ax)  (3.6 


Both  the  regular  and  singular  parts  are  antisymmetric  about 
Ax  -  0  .  The  regular  part  has  been  computed  and  is  shown  in 
Fig.  3.1*  It  has  been  checked  against  the  corresponding  function 
Wxq  (®xq  "  "  (R,g)  previously  tabulated  by  weissinger1^. 

As  in  the  case  of  (R^  ,  we  expand  (Rg  in  a  double  Fourier 
cosine  series  in  the  region  0  <  <t>  £  it  and  0  £  qr  £  ir  as 


p-0  q-0 


cos  p0  cos  q«p 


(3.7 


where  the  coefficients  J  ^  are  given  by  Eq.  (2.10)  with  (R^ 
replaced  by  Rg  and  h^q  by  Jp  q  •  Since  (Rg  is  an  odd 
function  in  Ax  ,  we  can  show  by  interchanging  4>  and  q>  that 


'p#q 


-  -  j 


q.p 


(3.8 


As  a  result,  Jq  0  '  ^1  1  '  ^2  2  '  all  identically  vanish. 

Also,  by  replacing  $  by  (tt-$)  and  q>  by  (*-9)  ,  we  find 


j  -  0 


if  (p+q)  -  even  integer 


(3.9 
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Using  the  relaxation  procedure  mentioned  earlier,  we  have 
computed  j  for  p,q  *  0  ,  1  ,  2  ,  ...7  for  X  =  1/4  , 

Xr  *  tL 

1/2  ,  3/4  and  1  until  the  difference  between  successive  compu¬ 
tations  was  less  than  +  1  x  10“'i  and  the  results  are  presented 
in  Tables  3*l-3-4.  Table  3*5  gives  a  comparison  between  the 

exact  value  of  &£  and  the  series  value  obtained  using  the  com- 

-4 

puted  jp  q  •  Agreement  is  to  within  +  5  x  10  ,  with  the 

series  value  for  the  smaller  X's  again  being  significantly 
better . 

With  the  expansion  of  into  a  singular  and  a  regular 

part,  the  direct-thickness  contribution  to  c  can  be  simplified 

J* 

to  matrix  form. 

3-3  Reduction  of  Direct-Thickness  contribution 

Since  each  of  the  two  factors  in  the  integrand  of  the  direct¬ 
thickness  contribution  is  composed  of  a  singular  and  a  regular 
part,  the  contribution  of  this  term  to  c^  could  be  reduced  to 
matrix  form  in  exactly  the  same  fashion  as  was  the  thickness- 
induced  camber.  That  is,  we  would  consider  the  product  of  the 
singular-singular,  singular-regular,  regular-singular  and  regular- 
regular  parts  and  evaluate  each  in  turn.  However,  it  is  more 
convenient  to  break  up  the  integration  in  the  following  manner. 

First  consider  the  product  of  t'  and  the  Cauchy  singular 
part  of  Zg  ,  or 


0.00795  0  -0.00001  0  0.00001  0  0.00001 
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0.01099  0.01139  0.01103 


o  r< 


O  X 


EXACT  VALUE  OF 
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(3.10) 


This  is  identified  as  the  thin-airfoil  theory  formula  for  the 
linearized  pressure  coefficient  of  a  two-dimensional  airfoil 
■action  with  thickness  t  and  no  camber,  see  for  example  Ref.  30. 
For  HACA  sections,  values  of  the  local  surface  velocities  have 
been  tabulated  in  Appendix  1  of  Ref.  36  and  it  is  better  to  use 
these  more  exact  values  multiplied  by  (-2)  rather  than  to 
evaluate  the  integral  of  Eq.  (3. 10). 

Since  has  been  expanded  in  the  same  form  as  ,  the 

remainder  of  the  integral,  or  the  product  of  t '  and  (Rg  ,  follows 
as  in  Sections  2.3  and  2.7  using  Jp  q  instead  of  h^  ^  ,  or 

X 

2  f  t'(i)  <R2(Ax)  df  -  -  {*)*[*]{«}  (3.11) 

-X 


Where  the  minus  sign  is  for  convenience  later  and 


(3.12) 


T  denotes  the  transposed  matrix,  and  [ T]  is  a  square  matrix  whose 
elements  are  functions  only  of  X  and  are  given  in  Table  3.6. 


2X 


(5  j0,l  +  35  J0,3  +  99  j0,5  +  195  J0,7 


("  Jl,0  +  15  j1.2  +  ft  jl.U  +  Ift  Jl.6  + 


,  1  .  1  .  l  4  1  ,  . 

5  J2,l  *  35  J2,3  "  99  J2.5  -  I95-  J2.7 

(-  j3.0  +  H  J3.2  +  ft  J3.4  +  I5J  J3.6  + 

I  J4.1  ‘  &  J4,3  ‘  ^  j4.5  ‘  135  J4.7  + 

(-  J5,0  +  i  J3.2  +  ft  J3.4  +  J5.6  + 

(-  J  j6.1  -  ^  J6.3  "  ^  J6.5  '  ik  J6.7  + 


••>/!  (-  Jl.O  +  1  jl,2  +  A-  J1.4  +  £  J1.6  +  •••) 


■  •oil 


■}li  (-  J3,0  +  5  J3,2  +  A"  J3.4  +  A  j3,6  + 


•ojl  (-  J5,o  +  5  J5.2  ♦  A  %,4  +  A  J5.6  +  •••) 

0/1 


("  I  jO.X  +  S'  jO, 


,1  ,  1  , 

(5  ^2,1  ■  S’  ^2, 


^3  j4,l 


(I  J6,l  " 


Him  Him 


^|ir\  Mlin  -tliA 


J0.3  +  A  j0,5  +  A  j0.7  + 


(-  5  Jl.C  -  A  j1.2  +  lA  J1.4  +  J1.6  +  •••)x2 


J2,3  "  A  J2,5  “  A  J2,7  +  •'■)X 


1  J3.C  “  A  J3.2  +  lA  J3.4  +  ?A  J3.6  + 


)* 


J4,3  '  21  J4.5  '  A  J4.7  + 


,)X 


5  J5.0  •  A  j5.2  +  eA  j5.4  +  lA  J5,6  +  *")x2 


j6.3  ‘  21  J6.5  "  A  J6.7  +  "'^X 


("  I  J0.1  +  AJ0.3  +  AJ0,5 


A  J2.1  '  A  J2,3  ”  A  j2,5 


A  J4.1  "  35  J4.3  ■  A  J4.5 


A  J6.1  -  A  J6,3  "  A  J6,5 


TABLE  3-6 
MATRIX  [ T ) 


5.3  +  h  J0.5  +  ft  J0,7  +  •")x3 


>.3  _  h  J2,5  ‘  ft  J2,7  +  •")x3 


(’  fJ0,l  -&Jo,3  +  ftJ 


0. 


5  J1.0 


-  35  Jl,2  +  ft  *1.4  +  Ilk  Jl.6  + 


(f  J2.1  +  J2,3  '  ft  J2. 


0  |  j3.0  “  ^  j3.2  +  ft  j3.4  +  IT5F  J3.6  +  •••)x4  0 

14,3  "  $3  ^4,5  "  ft-  J4,7  +  •••)x3  0  (7  J4, 1  +  Z?  *4,3  "  ft-  *4, 

0  ("  k  *5,0  ~  A  J5.2  +  ft  J5.4  +  lfk  J5.6  +  ••*)x4  0 


6,3  '  ^  J6,5  '  ft  J6.7  +  •”)x3 


0 


(f  *6.1  +  5J  j6,3  '  ft  J6, 


J0,l  -  55  j0,3  +  ft  J0.5  +  iSt  J0.i 


*7  J2,l  +  55  J2,3  ■  ft  j2.5  "  12-57  j2,l 


(7  J4,l  +  55  j4,3  *  ft  J4,5  ‘  I5&7  J4,1 


0 


(7  j6,l  +  is  J6.3  ‘  ft  j6,5  •  I5&T  j6.1 
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The  leading  forty-nine  elements  of  [t]  are  given  in  Tables 
3.7-3*10  to  an  accuracy  of  +  1  x  10  for  X  ■  1/4  ,  1/2  ,  3/4 
and  l  . 


3*4  Reduction  of  Effective-Camber  contribution 


As  noted  in  section  3-l»  the  integrand  of  the  effective- 
camber  contribution  is  already  in  a  form  suitable  for  integration. 
The  reduction  is  curled  out  by  decomposing  the  integrand  into  four 
parts  corresponding  to  the  singular -singular,  singular -regular, 
regular-singular  and  regular-regular  contributions  of  C  and  1^  . 
Bach  of  these  is  evaluated  in  turn  as  in  Sections  2. 4-2. 7  with  the 
aid  of  the  relationships  sin  nq>  •  sin  cp  -  ^{cos(n-l)<p  -  cos(n+l)<p) 
and  cot  Jjpp  »  (1+cos  <p)/sin  cp  and  the  following  definite  integrals 
«Aiich  are  derived  in  the  Appendix, 


7T 

J  in(coscp  -  cos$)2  dqp  -  -  it  in  4 


(3.13) 
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J  com  nqp  Jn(cosq>  -  cos#)2  dcp  ■  -  ~  com  n#  ,  n  >  1  ( 3 •  14 ) 


and 


7 r 
/ 


sin  vq>  cos  q#  sin  qp  dcp  -  g  ,  v  -  1  &  q  ■  0 


if  ,  v  -  q*l  4  q^l  (3.15) 


TABLE  3-8 


0.0000  0  -0.0001  0  0.0001  0  0.0000 


TABLE  3.10 
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Finally,  the  four  contributions  are  added  together  to  give  the 
Mtrix  fora  of  the  effective-caaber  contribution,  or 

tr 

k  J  c<<)^(®)  -  QliaC®))  -  {*)T[sHc}  (3.16) 

o 

Where  [ s]  is  given  in  Table  3.11.  The  leading  forty-nine 
elements  of  [s]  have  been  computed  to  an  accuracy  of  +  1  x  10""^ 
and  are  given  in  Tables  3.12-3.15  for  X  -  1/4  ,  1/2  ,  3/4  and  1  . 


3-5 


Final  Form  of  c 
_ _  _  £ 


The  final  form  of  the  duct  surface  pressure  coefficient  is 
now  obtained  by  substituting  Eqs.  (3.5) ,  (3*6),  (3>10),  (3. 11), 
and  (3.16)  into  Kq.  (3.3)  to  get 


cp  -  ±c(i)  +  F(x)  -  ($}*([T](a)  +  [site))  (3.17) 

Thus  we  have  succeeded  in  reducing  the  calculation  of  the 
aerodynamic  loading  on  an  annular  airfoil  to  a  simple  matrix 
operation  on  {a}  and  (c)  or  equivalently,  on  the  duct  sec¬ 
tional  properties  Aq  ,  ,  Ag  . . .  and  €q  *  €i  »  €2  '**  ainc® 

they  are  related  to  (a)  and  (c)  through  Eqs.  (2.6),  (1.20) 
and  (2.25).  The  matrix  method  of  computing  c  is  ideally  suited 

if 

for  engineering  calculations  and  preliminary  design  studies  since 
it  is  both  simple  and  rapid  and  the  effects  of  varying  any  or  all 


nX  (4  hQ  0  +  ^  tn 
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-n. 


-  4  h(1 
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(“  2  h4,0  +  h4,2^ 


2  h6,0  +  h6,2^ 
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h6,2  +  h6.4) 


TABLE  3.11 
MATRIX  [S] 


(-  h1>3  + 


(-  h3<3  +  h3(. 


h5.3  +  h5.: 
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0  (h^  2  "  ij)  0  (^o 1 4  ”  ^ 

hl,l  +  hl,3  _  &  °  (-hi  3  +  h15)  0  (-h1.5  + 

o  (- h2,2  +  h2,4  -  W5  0  (■  h2,4  +  h2,6)  0 

h3.1  +  h3.3  +  &  0  (-  h3  3  +  h3  5  -  £)  0  (-h3.5  + 

°  (-  h4,2  +  h4,4  +  0  (-  h4,4  +  h4,6  *  EF*  0 
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°  (“  h6,2  +  h6.4}  0  ("  h6,4  +  h6.6  +  ifc  0 


TABLE  3.11 
MATRIX  [S] 
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of  the  parameters  t  ,  e  ,  and  X  may  be  easily  determined.  An 
estimate  of  the  error  in  any  calculation  may  be  made  from  a 
"reasonable  guess"  at  the  magnitude  of  the  elements  Which  have 
been  neglected  in  the  tabulated  matrices.  Generally,  the  error 
will  increase  as  X  increases. 

3-6  Sectional  Radial  Force  and  Pitching  Moment 

For  an  annular  airfoil  at  zero  incidence  in  a  perfect  fluid 
the  net  forces  and  moments  acting  on  the  airfoil  vanish.  How¬ 
ever,  if  we  consider  a  section  of  the  duct,  then  we  find  that 
both  a  radial  force,  or  lift,  and  a  pitching  moment  are  present. 
The  form  of  the  corresponding  force  and  moment  coefficients,  cf 
and  cm  respectively,  are  the  same  as  for  a  two-dimensional 
airfoil,  or 

-  *(c0  +  hPx)  (3-18) 


§■  (cx  -  c2)  (3.19) 

radial  force  is  positive  outwards  and  the  pitching 
taken  about  the  quarter  chord,  positive  nose-down. 


and 


cm  - 


where  the 


CHAPTER  POUR 


NUMERICAL  EXAMPLES 


4.1  Illustrative  Calculation 

In  order  to  illustrate  the  procedure  for  computing  the 
duct  pressure  distribution  we  will  work  out  an  example  in  detail, 
we  take  an  NACA  0008  airfoil  with 


c  ■  -  tan  5* 

X  -  1/2  (4.1) 


First,  we  determine  (e)  .  It  follows  immediately  by  writing  c 
in  the  form  of  Eq.  (2.1)  and  is 


{«) 


r  > 

-0.1750 

o 

<  o  > 


(4.2) 


Next,  for  the  NACA  0008  airfoil  5(tmax/c)  -  0.4  ,  and  from 
Table  2.1  we  find  for  \  ■  1/2  , 
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(a)  ■ 


0.1188 

-0.2521 

-0.1238 

0.1951 

-0.3248 


(4.3) 


To  obtain  (€fc)  we  multiply  (a)  by  the  appropriate  matrix 
[q]  from  Table  2.12  to  get 


0.0005 

-0.0356 

-0.0116 

0.0009 

-0.0005 

0.0002 

-0.0002 


(4.4) 


for  the  first  seven  elements . 


The  two-dimensional  Glauert  coefficients,  or  (c)  -  (c^)  ,  are 


(c2“D) 


-0.1755 

0.0356 

0.0116 

-0.0009 

0.0005 

-0.0002 

0.0002 


(4.5) 


I 
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and  so  the  Glauert  coefficients  of  C  are  obtained  from 
Sq.  (1.20)  and  Table  1.2  a> 


Cc} 


-0.1890 
0.0088 
0.0070 
<  -0.0009  > 
0.0005 
-0.0002 
0.0002 


(4.6) 


Note  that  the  three-dimensional  effect  of  the  duct  changes  only 
the  first  three  two-dimensional  Glauert  coefficients. 

With  the  required  Glauert  coefficients  determined,  the 
expression  for  the  surface  pressure  coefficient,  Eq.  (3. 17),  may 
be  simplified.  The  first  term,  +_  C  ,  is  found  from  Eqs.  (1.19) 
and  (4.6).  The  contribution  of  F  is  obtained  from  the  basic 
thickness  form  data  in  Appendix  1  of  Ref.  26.  The  matrix  opera¬ 
tions  [T]{a)  and  [s](c)  are  carried  out  using  Tables  3.8  and 
3.13  and  Eqs.  (4.3)  and  (4.6),  yielding 


tT](a)  +  (s]{c) 


-0.2195 

0.0945 

0.0036 

<  0.0007  > 

0.0000 

0.0000 

0.0000 
V-  J 


(4.7) 
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With  this  information  and  (0)  from  Eq.  (3*12),  cp  has  been 
calculated  and  is  plotted  in  Fig.  4.1.  For  comparison  purposes, 
the  effect  of  doubling  X  While  keeping  t  and  e  the  same 
has  also  been  computed  and  is  shown  in  the  same  figure. 

The  effect  on  c^  of  varying  t  is  given  in  Fig.  4.2  for 
c  ■  -  tan  5°  and  X  -  1/2  .  comparison  of  the  results  with  the 
previous  figure  reveals  that,  for  the  particular  configuration 
investigated,  doubling  the  chord-to-diameter  ratio  has  a  much 
smaller  effect  on  cp  than  a  nearly  corresponding  increase  in 
the  thickness  ratio. 

4.2  Comparison  With  Previous  Results 

To  assess  the  accuracy  of  this  new  method,  we  have  calcu¬ 
lated  the  velocity  and  pressure  distributions  on  several  annular 
airfoils  to  compare  with  the  experimental  and  theoretical  results 
of  previous  investigators.  The  only  wind  tunnel  measurements  of 
pressure,  or  velocity,  distributions  published  to  date  tfiich  we 
have  found  are  those  of  the  Bureau  Technique  Zborowski11.  in 
their  tests,  the  circumferential  variation  of  the  inner  and  outer 
surface  pressures  on  an  NACA  66-006,  X  ■  1  ring  airfoil  was 
measured  at  several  angles  of  attack.  Their  data  which  is  appro¬ 
priate  to  the  present  study,  i.  e.,  the  zero  angle  of  attack 
results,  has  been  replotted  to  within  ±  0.005  in  Fig.  4.3  and 
compared  with  the  distribution  as  calculated  from  Eq.  (3 .17)- 


FIGURE  4.1 


VARIATION  OF  PRESSURE  DISTRIBUTION 
WITH  DUCT  CHORD-TO -DIAMETER  RATIO 


Except  in  the  immediate  vicinity  of  the  leading  and  trailing 
edges,  the  results  differ  by  less  than  10 #.  Better  agreement 
is  obtained  for  the  inner  surface  distribution. 

The  pressure  distribution  on  a  conical  ring  airfoil  of  zero 
thickness  as  computed  from  the  rheoelectric  tank  analogy  results 
of  Malavard12  is  compared  with  that  of  the  present  theory  in 
Pig.  4.4.  Since  X  =  2/3  for  this  configuration,  we  have  used 
values  of  the  necessary  matrices  which  were  estimated  by  linear 
interpolation  of  the  results  for  X  =  1/2  and  X  =  3/4  •  Satis¬ 
factory  accuracy  can  be  achieved  in  this  manner.  In  this  example, 
good  agreement  is  obtained  for  the  outer  surface  distribution  only. 
The  discrepancy  between  the  inner  surface  results  has  not  been 
resolved . 

Bagley,  Kirby  and  Marcer^  have  developed  a  method  for  finding 

the  velocity  distribution  over  annular  airfoils  in  incompressible 

flow.  The  basis  of  their  work  is  the  approximation  of  certain 

31 

integrals  by  Weber  "sum -functions"  of  the  airfoil  ordinates  and 
the  satisfaction  of  the  resulting  set  of  n  linear  equations  at 
n  points  on  the  airfoil  chordline,  comparison  with  their  results 
for  the  thickness -induced  radial  velocity  (less  the  term  +  f) 
on  an  RAE  101-10#  ,  X  =  1/2  ring  airfoil  is  shown  in  Fig.  4.5* 
Excellent  agreement  is  obtained,  the  results  being  indistinguishable 
except  near  the  trailing  edge.  Fig.  4.6  compares  the  induced  axial 
velocities  on  an  RAE  101-10#,  X  -  1  ring  airfoil;  again,  rela¬ 
tively  good  agreement  is  obtained. 
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We  have  found  that  the  computation  of  the  pressure  distri¬ 
bution  at  thirteen  chordwise  stations  takes  approximately  one 
hour  on  a  desk  calculator  once  (cj  and  (a}  are  determined. 
The  advance  tabulation  of  {$)  corresponding  to  these  stations 
proves  to  be  of  great  assistance  in  facilitating  the  numerical 
work. 


CONCLUSIONS 


A  simple,  rapid  method  for  the  determination  of  the 
linearized  surface  pressure  distribution  on  streamlined  ducted 
bodies  In  axial  flow  has  been  developed.  The  pressure  coeffi¬ 
cient  Is  composed  of  camber  and  thickness  terms  analogous  to 
two-dimensional  flow  plus  an  additional  term  expressed  as  a 
matrix  operation  on  the  duct  sectional  properties.  Tables  of 
the  necessary  matrices  for  four  representative  duct  chord-to- 
dlameter  ratios  are  given;  matrices  obtained  by  Interpolation 
from  these  values  can  be  used  for  "in-between”  chord -to-diameter 
ratios . 

The  method  is  ideally  suited  for  engineering  calculations 
and  preliminary  design  studies  and  the  effects  of  varying  any 
or  all  of  the  thickness  and  camber  distributions  and  chord-to- 
diaraeter  ratio  may  easily  be  determined. 

Comparison  of  predicted  pressure  and  velocity  distributions 
for  several  annular  airfoils  gives  good  agreement  with  previous 
experimental  and  theoretical  results.  A  typical  numerical  com¬ 
putation  for  a  given  configuration  may  be  done  on  a  desk  calcu¬ 
lator  in  about  one  hour. 

The  theory  may  be  extended  to  include  the  interference 
effects  of  a  centerbody  and/or  a  propeller  as  well  as  duct  angle 
of  attack. 
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APPENDIX 


Derivations  o£  certain  integrals  referred  to  in  the  main 
portion  of  the  text  are  given  here* 

it 

cos  n <t>  Jn(l+cos0)  d <t>  ,  n  >  1 

0 

He  integrate  by  parts  with  u  =  in(l+cos<£)  and 
dv  ■  cos  n<t>  d <t>  .  Upon  carrying  this  out,  we  find  that  the  con¬ 
tribution  of  uv  vanishes  and  we  are  left  with 

it 

1  f  sin <t>  sin  n$ 

n  J  -  d* 

0 

This  is  evaluated  in  turn  for  each  n  by  expanding  sin  n<J>  in 
powers  of  sin  <t>  ,  removing  the  factor  sin  <t>  ,  and  dividing 
(l+coa$)  into  sin20  -  (l+cos</>)  (l-costf>)  .  The  integral  is  then 
readily  reduced  and  the  results  generalized  to  give 

J  cos  n0  in(l+cos0)  d 0  ■  - —  t  n  >  1  (A.l 

0 
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/•***£»* 

0 

we  integrate  by  parts  with  u  «*  In  and 

dv  ■  sin*  d0  .  The  contribution  uv  vanishes  and,  after 
trigonometric  reduction,  v  du  is  simply  2  ,  so  that 

/  •**  *  -*  -  *■  (x-2) 
0 


ir 

J'  cos2n*  sin*  In  d*  ,  n  >  1 

0  _ _ 

We  integrate  by  parts  with  u  -  In  ®"d 

dv  «  cos2"*  sin*  d*  .  Then  uv  vanishes  and  we  are  left  with 


2 

5n+T 


J'  COS2"*  d* 


Changing  to  a  period  of  tt/2  ,  we  can  integrate  this  in  turn  by 
means  of  Formula  858*44  of  Ref.  20  to  obtain 


f  cos 2n*  Sin*  In  d*  - 


(2n-D 

,  2n 


n  >  1 


(A.3) 


71 

J  ln(coscp  -  cos#)2  d«p 


We  integrate  by  part*  with  u  -  ln(coscp  -  cos# )  and 
dv  -  dtp  to  get 


71 

w  ln(l+cos#)2  +  2  J' 


coscp  -  cos? 


This  integral  is  somewhat  similar  in  form  to  the  Glauert 
integral  and  can  be  evaluated  by  complex  contour  integration. 
Making  the  substitution  z  ■  e^  ,  we  study  a  related  integral 
of  the  form 


%  (*-e  *)(z-e 


(z-eiV)(z-e-iV) 

where  c  is  the  closed  curve  composed  of  ci  '  c2  '  C3  '  the 
indentations  at  +  #  ,  ,  and  Lg  ;  viz . , 
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The  integrand  Is  analytic  inside  the  region  bounded  by  c  and 
thus,  using  the  residue  theorem^2  of  complex  variables,  we  have 

/  ♦  /_* +  /  +  /  ♦  /  ♦  /  +  /  * 0 
C1  e  C2  ®  c3  L1  C4  L2 

On  the  identation  at  +  0  ,  z  ■  e^  +  ceia  Where  a  turns 

through  an  angle  of  -  n  .  The  contribution  here  is  evaluated 

taking  the  limit  as  e  ■*  0  .  At  the  lower  indentation,  -  0  , 

_i<i  in 

z  =  e  +  e  ,  and  after  passing  to  the  limit  €  ■+  0  ,  we  find 
that  the  contribution  here  exactly  cancels  that  at  +  0  .  On 
,  z  *  re  while  on  Lg  ,  z  =*  re  .  The  contributions  of 
,  and  Lg  are  combined  and  simplified  to  yield 


f  .r  r  -tr-*:1) 

*L  r**  2r  cos0  +  1 


it  in  4(l+cos0)  +  2ir  in  6 


On  ,  z  ■  deitf  ,  -  *  6  £  ir  ,  and,  after  simplification,  we 

find  that 


-  2ir  Jn  4 


collecting  these  results,  we  have 


■  -  ir  in  4(1+00*0)* 


Cj+Cg+Cg 
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■in  vq  cos  q<p  s inq  dq 

0 

We  decompose  the  integrand  into  the  sum  of  four  cosines  by 
Formula  402.05  of  Ref.  20,  or 

4  sin  vq  cos  qq  sinq  -  [cos(q+l-v)q  +  cos(q+l+v  )q 

-  cos(q-l+v)q  -  cos(q-i-v  )<p] 

me  subsequent  integrations  are  easily  carried  out,  giving 

IT 

f  sin  vq  cos  qf  sinq  dq  -  £  ,  »  ■  1  fc  q  ■  0 

0 

-  +  J  #  v-q+lftq^l  (A. 7) 


